Decohering localized waves 
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In the absence of confinement localization of waves takes place due to randomness or nonlinear- 
ity and relies on their phase coherence. We quantitatively probe the sensitivity of localized wave 
packets to random phase fluctuations and confirm the necessity of phase coherence for localization. 
Decoherence resulting from a dynamical random environment leads to diffusive spreading and de- 
stroys linear and nonlinear localization. We find that maximal spreading is achieved for optimal 
phase fluctuation characteristics which is a consequence of the competition between diffusion due 
to decoherence and ballistic transport within the mean free path distance. 
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I. INTRODUCTION 

Ever since the absence of diffusion of a quantum par- 
ticle due to randomness on a lattice was predicted by 
P.W. Anderson wave localization has become an in- 
tensively studied phenomenon. Of particular interest is 
the wave character of the quantum objects which localize. 
In contrast to particles waves are fields characterized by 
amplitude and phase, and can tunnel under a potential 
barrier but also be back-scattered above a given barrier. 
In the absence of confinement wave localization implies 
that the wave does not escape to infinity even if it was 
energetically allowed. Consequently Anderson localiza- 
tion is an interference phenomenon which relies on phase 
coherence. 

Localization on a lattice also occurs as result of an ap- 
plied DC field where it leads to Bloch oscillations 0, y| . 
It can even show up in translationally invariant lattices 
in the form of discrete breathers as a result of nonlinear 
interactions [111. Another source for wave localization 
can be a quasi-periodic potential interpolating between 
uncorrelatcd disorder and perfect periodicity. In this case 
localization has been predicted by Aubry and Andre [6[ 
only when the strength of the quasi-periodic potential ex- 
ceeds a certain critical value. The latter case is largely de- 
bated. Initially the additional importance of high enough 
incommensurability was stressed 04^1 • More recently, 
the occurrence of a metal-insulator transition has even 
led to the claim that localization in the Aubry-Andrc 
model has particle character [lfj , and therefore the phase 
coherence should not matter. 

In this work we probe linear and nonlinear localiza- 
tion through its sensitivity to phase coherence. We find 
that loss of phase coherence directly relates to delo- 
calization. In experiments decoherence generally arises 
due to random temporary fluctuations as a result of in- 
evitable coupling of the ideal system to its environment 
[ill [l2j . Therefore our considerations are also related to 
the conceptual understanding of the effects of decoher- 



ence which is assumed to play a key role in the translation 
of the quantum world to the "classical" pict ure obtained 
through generic measurement processes [13|. 

The importance of minimizing fluctuations becomes 
evident from the decade-long difficulties of experimen- 
tally observing Anderson localization. In disordered elec- 
tronic systems, where it was initially predicted, localiza- 
tion is hindered by phonon and electron assisted variable 
range hopping [1J, [l5|]. Anderson localization has also 
been related to the study of quantum chaos in the quan- 
tum kicked rotor system [ll|. Here, classical chaotic dif- 
fusion is suppressed by quantum interferences and leads 
to localization in momentum space. In this system the 
destructive impact of noise and decoherence on Anderson 
localization has been discussed as well, both theoretically 
[ItJ and experimentally with ultracold atoms in a pulsed 
standing light wave [HI, [l9j]. However, direct observa- 
tion of Anderson localization of matter waves has been 
achieved only recently [20|, [2l[ with the creation of Bose- 
Einstein condensates where fluctuations are sufficiently 
reduced. Apart from quantum systems, the wave char- 
acter of Anderson localization, with the inherent role of 
coherence, has also been demonstrated with light waves 
in coupled optical wave guides [23 |. 

In the following we will investigate wave packets which 
are localized not only as a consequence of random dis- 
order, but also of quasi-periodic disorder and DC fields 
(linear localization) and of interactions (nonlinear local- 
ization). We will show that specifically fluctuations in 
the phases, resulting from coupling to a random environ- 
ment, lead to delocalization which is directly related to 
the loss of coherence. The loss of localization is observed 
as a generic onset of a diffusive spreading regime. Inter- 
estingly, maximal spreading and therefore delocalization 
is not achieved for the strongest or most frequent dephas- 
ing. There exists an optimal rate and strength of dephas- 
ing which maximizes the wave packet's extent at a given 
time. The optimal rate can be derived from knowledge 
of the localization volume of a given wave packet. All re- 



2 



suits also remain true for the case of quasi-periodic dis- 
order, therefore localization in the Aubry-Andre model 
has, as in the Anderson localization problem, pure wave 
character and in the case of matter waves is of quantum 
origin. 



II. LOCALIZED WAVE MODELS 

We consider the model of the one-dimensional discrete 
nonlinear Schrodinger equation 

d 

i -Q t t Pi = e i^i + V'J+i + ipl-x + P\~>Pi\ ipi, (i) 

where ipi is a complex field at site I with on-site energy ei 
and nonlinearity strength (3. It can be derived from the 
Hamiltonian 

H = J2 *l\M 2 + + W+i^l) + fl^l 4 (2) 

i z 

by 0; = dH/d(iipf). Varying the total norm S := 
I "0(1 2 is equivalent to varying /3, therefore it is always 
possible to fix S = 1. Then \ipi\ 2 can be identified as 
the norm density on site I. Both the energy H and the 
norm S are integrals of motion. The discrete nonlinear 
Schrodinger equation allows for the investigation of all 
localization phenomena described above. 

Linear localization of wave packets is obtained for (3 = 
and suitable potentials e/ ^ 0. The resulting set of 
linear differential equations may be decoupled in terms 
of normal modes (NMs). Separation of variables {tpi(t) = 
Aie~ lXt ) leads to the eigenvalue problem 

XAi = eiAi + Ai+i + Ai-i (3) 

for the only site dependent (in general complex) ampli- 
tudes Ai. The normalized eigenvectors A v ^ are the 
NMs having cigcnfrequencies A„. For localized NMs any 
initial excitation containing only a finite number of them 
will stay localized in time, leaving the NMs oscillating 
with their eigenfrequency respectively but independently 
from each other. As will be shown later, real amplitudes 
A/,z (especially including negative ones) ensure localiza- 
tion of a single NM. 

We consider the following three cases of on-site ener- 
gies which yield localized NMs. (I) Disorder (Anderson 
model) : the on-site energies are random and chosen uni- 
formly from the interval [— where W denotes the 
disorder strength. The eigenfrequencies A„ lie in the in- 
terval [-2 - ^,2 + ^]. The NMs decay exponentially 
as A u j ~ e~~ 1 ^ with a localization length £(A„). (II) 
DC Field (Wannier-Stark ladder): the field E determines 
the linear growth of the on-site energies, e/ = El. The 
eigenfrequencies are X v = Ev and the NMs are given by 
Bessel functions of first kind Ji(x) as A u j = Ji_ v (2/E) 
(Ill) Quasiperiodic potential (Aubry-Andre model): 



the on-site energies satisfy e; = £cos(27ra/). The com- 
mensurability is characterized by the irrational parame- 
ter a while the parameter £ describes a relative strength, 
similar to W in the Anderson model. A condition for 
localized NMs is that a be as far as possible from a ra- 
tional number , where a standard choice is to take the 
inverse golden mean a — (23j . Then NMs are ex- 

ponentially localized for ( > 2 in real space and for ( < 2 
in Fourier space (self-duality) . Consequently there exists 
a metal-insulator transition at £ = 2 in either space. 

Adding a nonlinearity in general destroys the integra- 
bility of the linear system by inducing frequency shifts 
which result in excitations of overlapping NMs [24l HH . 
For localized NMs small enough nonlinearities have been 
shown to cause derealization with subdiffusive spread- 
ing [24|, [111 [27| . In contrast, a large enough nonlinear- 
ity can lead to localization in form of discrete breathers 
when the nonlinear frequency shifts exceed the linear 
eigenfrequency spectrum (self-trapping) d, [28| . Discrete 
breathers are exponentially localized and their time- 
dependence is characterized by a single frequency Of, as 
i)i{t) = Aie 1 ^. 

To study nonlinear localization independent of linear 
localization, we consider Eq. ([IJ with e; = 0,V7,/3 > 0. 
Then the NMs of the underlying linear system (Eq. Q 
with ei = 0) are extended and localization can only be at- 
tributed to discrete breathers. Any not self-trapped part 
of a wave packet will spread ballistically over the entire 
chain. Note that overlap of several discrete breathers in 
general does not lead to new localized wave packets with 
a quasiperiodic time-dependence @. 

To quantify the spatial extent (localization volume) of 
an arbitrary wave packet we compute the participation 
number P = \Ai\ 4 which measures the number of 

considerably excited sites [2!|, and the second moment 
(variance) m,2 = X)i( m i — 2 |^z| 2 which measures the 
average norm spread around its center mi = J2i l\Ai\ 2 . 
Since mi has the units of a square distance, a localization 
volume in one dimension is defined proportional to yJrrUi 



III. LOCALIZATION AS A COHERENT 
EFFECT 

In this paper we want to establish a direct relationship 
between coherence and localization. The main goal is to 
show how temporal phase fluctuations always destroy lo- 
calization. Consequently phase coherence is necessary for 
localization, which is a manifestation of its crucial wave 
character. Considering localization in form of localized 
NMs and discrete breathers, we will first identify max- 
imally localized wave packets with maximal coherence. 
Then we will investigate the effects of phase fluctuations 
and decoherence. 

To characterize a wave packets' degree of coherence we 
will utilize the basic property of the second order com- 
plex degree of coherence 7 [3l[ which is closely linked to 
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the visibility of interference fringes in related interference 
experiments. The essence of second order coherence lies 
in the investigation of correlations between two space- 
time points. It is furthermore a prerequisite for coher- 
ence of higher orders that contain correlations of more 
space-time points. The second order complex degree of 
coherence is defined as [3l| 

y(h h t) •= r (^,r) f 4 ) 

which is a normalized version of the mutual coherence 
function 

T 

T(h,l 2 ,T) = lim i ( ri(t)Mt + r)dt. (5) 

T->oo 1 J 


Complete coherence is achieved for |7(ii, h, t)\ — 1, while 
partial coherence corresponds to < \j{h, h, t)| < 1 
and complete incoherence to (7(^1, 121 t)\ = 0. Note that 
the averaging in Eq. (|S|) assumes statistically stationary 
evolving fields ipi,ip2- 

Inserting a generic solution of the form Aie lujt with site- 
independent lu and time- independent amplitudes^; into 
Eq. ((3]), it is straightforward to obtain I7I = 1. Thus a 
single NM or discrete breather always possesses complete 
coherence, when uj is identified with either X u or flf,, re- 
spectively. The phases at different sites of these coherent 
solutions remain locked together in time. It can also be 
shown easily that a superposition of more than one NM 
will only be partially coherent. Since nonlinear localiza- 
tion is only expected in form of discrete breathers it is 
in general always fully coherent. Note that here we have 
not made any restrictions to the complex amplitudes Ai 
except for time-independence. Especially we have not re- 
quested that the solutions (e.g. the NMs) be localized. 
Therefore coherence is not a sufficient condition for lo- 
calization. 

In order to obtain a precise necessary condition be- 
tween localization and coherence let us consider the evo- 
lution of the norm density pi — ipt^l — \Ai\ 2 - Norm 
conservation of Eq. ([T]) ensures that the continuity equa- 
tion 

o=^ + (im-i0. (6) 

holds. Inserting Eq. (TTJ) into Eq. © one obtains the 
familiar result that a phase difference in a complex field 
ipl = \^i\e- lipi leads to a norm density current 

3i = -2\i>l+l\\i>l\ sin(^z+i - <pi). (7) 

For time-independent pi the first term in Eq. (jSJ) has 
to vanish, requiring ji to be actually independent of site 
I, Now consider that localization implies \Ai\ — > for 
I —> ±00. Then the sites with \Ai\ =0 demand ji = V7 
in Eq. (JTJ). Thus for any localized solution ji need not 
only be constant in space but even strictly equal to zero. 



The sites where \Ai\ ^ can only satisfy ji = when 
(pi+x — <pi= TO7r, to = 0, ±1, ±2, ... Consequently the con- 
dition of locking together the phases of neighboring sites 
at differences of a multiple integer of it is necessary for 
localization of a single NM and discrete breather. This 
is equivalent to choosing real amplitudes Ai and espe- 
cially allowing negative values as well. Since the most 
localized wave packets consist of either a single NM or 
discrete breather in a linear or nonlinear model respec- 
tively the above condition is also necessary for maximal 
localization. 



IV. DECOHERING AND DELOCALIZING 

Now we will quantitatively analyze the impact of tem- 
poral phase fluctuations which in general occur when cou- 
pling to a stochastic environment is considered. It is clear 
that for a maximally localized wave packet this will de- 
stroy the fixed phase relation of neighboring sites and 
lead to a local norm density current. Then complete co- 
herence is lost and spreading takes place. However, the 
wave packet does not have to delocalize completely as 
long as some partial coherence remains. We will numeri- 
cally show that an initial dephasing induced loss of coher- 
ence directly relates to the loss of localization and that 
persistent fluctuations finally delocalize the wave packet 
completely. 

We integrate Eq. JT]) using a symplectic SABAi inte- 
grator described in [321 ] . A dynamical random environ- 
ment is included by considering dynamical on-site ener- 
gies e/—>e/ + £;(£) where e\(t) is a random process. It 
can be easily shown that £/(i) defines a dephasing term 
in the equations of motion of the phases ifi . Numerically 
the dephasing is implemented by altering the phases be- 
tween integration steps of the unperturbed equations of 
motion. 

We choose two different dephasing schemes. The first 
one is a complete random dephasing which allows to 
observe the effects of dephasing as clearly as possible. 
Therefore the phases <fi are replaced by completely new 
random phases chosen uniformly from [0, 2tt] at certain 
times of the integration. Between these kicks on the 
phases the dynamics of a wave packet is solely governed 
by the linear or nonlinear equations of motion. Each 
kick can be considered as defining a new initial configu- 
ration for the integration where all the information of the 
old phases is lost. The second scheme is a quasiperiodic 
dephasing, similar as in [l2j . Its purpose is to probe the 
sensitivity of localization to only small phase fluctuations 
which are even not completely uncorrelated. The new 
phases result from the old ones as fi(t) — ipi + 6sin(/i;t) 
with a time-independent frequency pn chosen on each site 
randomly and uniformly from [0, max{pi}] and a con- 
stant strength of dephasing b. The frequencies pi are in 
general incommensurate but fixed during the integration. 
This leads to uncorrelated fluctuations between different 
sites and temporarily correlated fluctuations on each site. 
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FIG. 1: Comparison of (a) |7| and (b) 7712 in the Anderson model after 
complete random dephasing of a NM. Dcphasing is switched off after a 
number k of kicks: k— (black line), 1 (red), 2 (green), 3 (blue). In (b) 
only times arc shown when the wave packet is on average not spreading 
anymore. A direct correspondence between decohcrence and dereal- 
ization is observed when changes in the averages of |7(t)| (dashed lines 
in (a)) and Tnzft) (dashed lines in (b)) are compared. 



FIG. 3: Diffusive spreading in linear chains (Anderson model) for per- 
sistent dcphasing of initially one normal mode, (a) Complete random 
dcphasing with a time A between phase kicks of A— 10 red (r), 1 green 
(g), 0.01 blue (b) and 0.001 brown (br). The black (bl) line shows the 
evolution without dcphasing for reference, (b) Quasipcriodic dcphasing 
with strength b = 0.0001 black (bl), 0.0006 red (r), 0.01 green (g), 0.06 
blue (b) 1.1 brown (br). 




(a) number of kicks (b) number of kicks (c) number of kicks 



FIG. 2: Time averages of I7I (lower panel), 1712 (upper panel) in the 
(a) Anderson model, (b) Wannicr-Stark ladder, (c) Aubry-Andrc model 
when dephasing is switched off after a small number k of kicks. The 
averages for k— 1, 2 and 3 kicks in (a) correspond to the dashed lines in 
FigE Up to a number of 8 kicks (in (a), (c), 9 kicks in (b)) an increase 
(decrease) in 7 corresponds to a decrease (increase) in 1712. 



Here, the phase kicks are performed after each step of in- 
tegration. To calculate the complex degree of coherence 
we switch off dephasing after a certain number of kicks 
on the phases in order to obtain statistically stationary 
evolving fields. We use the time evolution at the center 
of norm and one neighboring site. 

Dephasing of a NM (FigO} shows that the decoherence 
is in general non-monotonic with respect to the growing 
number of phase kicks, however, a decrease (increase) in 
coherence similarly leads to a decrease (increase) in local- 
ization. The reason is a selective excitation and damping 
of overlapping NMs. Therefore the direct correspondence 
between decoherence and delocalization can be found in 
all linear models with localized NMs (FigJSJ). It holds 
for a small number of phase kicks as long as the excited 
NMs contain the two sites which are considered in the 
calculation of | — ^ | . 

In contrast, after a few phase kicks on a discrete 
breather, cither only a small part of the norm is radi- 
ated and the rest can again self-trap to form a new fully 
coherent wave packet, or the discrete breather delocal- 
izes completely. However, the approach to a new time- 
periodic trajectory of a discrete breather in phase space 
can take up arbitrarily long times Q. In this regime we 
observed seemingly localized (over accessible integration 



times) wave packets with only small deviations from com- 
plete coherence. The deviations seem to be larger when 
the initially radiated norm (and therefore delocalization) 
has been larger, approximately confirming a correspon- 
dence between decoherence and delocalization. 

A common result for both linear and nonlinear mod- 
els is that coherence on average decreases for a large 
growing number of phase kicks. This situation is close 
to an experimental one where phase fluctuations exist 
persistently. We observe that decoherence by persis- 
tent dephasing leads to complete delocalization and the 
generic onset of a diffusive spreading regime, 7712 ~ t and 
P ~ Vi (FigEland Figg]), as can be expected for sim- 
ilar dynamically disordered tight-binding Hamiltonians 
[33j . Note that the loss of norm from a discrete breather 
is properly characterized by the participation number. 
Here, the diffusive regime is preceded by jump-like in- 
creases of P which correspond to radiation of substantial 
parts of norm as small amplitude waves. It is expected 
that discrete breathers remain robust upon radiation of 
small amplitude waves [34l ] which can be treated as linear 
background of the high amplitude breather [35[ . Conse- 
quently, when dephasing is only performed very rarely, 
new self-trapping may occur leading to a series of step- 
wise increases of P (FigHJa)) . 

In general, the transition times before the asymptotic 
diffusive regime is assumed can vary. Especially for the 
quasiperiodic dephasing it may be expected that even 
for very small b normal diffusion commences after suf- 
ficiently long waiting time (see for the Anderson 
model). This shows the highly sensitive dependence of 
localization on phase coherence and minimal phase fluc- 
tuations. Moreover the results do not depend on the spe- 
cific range [0,max/x;] of dephasing frequencies fii. Con- 
sequently, in experiments localization is at best an inter- 
mediate regime, since arbitrarily small fluctuations can 
destroy localization. 

Another interesting effect of random dephasing is that 
increasing the rate or strength of dephasing does not nec- 
essarily lead to a decrease in localization and to faster 
spreading. It is rather observed that very strong and fre- 
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FIG. 4: Diffusive spreading for persistent dephasing of a discrete 
breather. The crossover from a ballistic to diffusive spreading regime 
is seen as a jump-like increase of P. (a) Complete random dephasing 
with time A between phase kicks of A— 900 red (r), 500 green (g), 
10 blue (b), 0.1 brown (br). The black (bl) line shows the evolution 
without dephasing for reference, (b) Quasipcriodic dephasing with 
strength b = 0.0001 black (bl), 0.00015 red (r), 0.0002 green (g), 
0.00025 blue (b), 0.0006 brown (br). 




0,5 1 1,5 2 
time between kicks 



0,003 0,006 0,009 0,012 
strength of dephasing b 



FIG. 5: Dependence of the diffusion constant D on (a) rate and 
(b) strength of dephasing in the Anderson model (black lines (bl)), 
Wannicr-Stark ladder (red lines (r)) and Aubry- Andre model (green 
lines (g)) with the asymptotic dependence in the insets. The occur- 
rence of a global maximum in each model is observed. Each data point 
of D results from linear regression of an average of rri2(t) over 100 real- 
isations of sequences of random dephasing. The connecting lines guide 
the eye. 



quent dephasing even may suppress the onset of normal 
diffusion for an increasingly long time. As a result, the 
wave packet extent at the final time of integration be- 
comes maximal for a certain optimal rate and strength 
of dephasing. 

When the maximal extent corresponds to a diffusive 
spreading regime it can be attributed a maximal diffusion 
constant D. This constant then also defines the optimal 
dephasing parameters for maximal derealization at later 
times. Note that in the nonlinear model (Fig0]) this is not 
the case. A discrete breather is destroyed best when the 
onset of normal diffusion can be avoided for a long time. 
Therefore optimal rates and strengths are different for 
different final times considered and cannot be identified 
with a macroscopic constant. 

For the linear models the diffusion constants are shown 
in Fig. [5j Apart from some local maxima there clearly 
exists a global maximum (see also [36j, |37[ for the An- 
derson model) in a wide range of rates and strengths 
of dephasing. Its nature can be revealed with the fol- 
lowing estimate in the Anderson model. Consider that 
the localization l eng th corresponds to a mean free path 
of quasiparticles [38| . An optimal rate would leave the 



wave packet spread ballistically over the new accessible 
localization length £ before the next phase kick is applied. 
Since after one kick only those NMs can be excited con- 
siderably, whose center of norm lies inside the localization 
volume (determined by £) of the initial wave packet, at 
most one localization volume becomes additionally ac- 
cessible. An upper estimate of the localization volume 
of a NM is of the order of 330/W 2 corresponding to a 
localization length of 100/ VF 2 [3{|. For W = 5 this gives 
a mean free path for quasiparticles of ~ 4. With a max- 
imal particle velocity (group velocity in Fourier basis) of 
two sites per time unit an optimal rate corresponds to 
one kick every two time units. This upper estimate is 
in good agreement with the numerically obtained value 
(Fig. Ufa)) of one kick per 1.6 time units. As a result, 
one could actually expect a ballistic spreading when the 
kicks match the times for completing the ballistic spread- 
ing into a new localization volume. However, we always 
observe the onset of a diffusive regime since the new lo- 
calization volumes variate while the rate of dephasing is 
kept constant. In consequence of the above, optimal pa- 
rameters of dephasing have to balance on the one hand 
decoherence and derealization of the initial wave packet 
and on the other hand should not suppress too much 
a possibly super-diffusive spreading by imposing normal 
diffusion. Therefore it becomes also clear that the opti- 
mal dephasing of a discrete breather cannot correspond 
to a diffusive regime since stepwise radiation of small am- 
plitude waves allows for ballistic transport to arbitrary 
distances in this case. 



V. CONCLUSIONS 

In conclusion, we investigated the effects of decoher- 
ence on linear and nonlinear localized wave packets and 
showed that phase coherence is a necessary condition 
for localization. Therefore localization, especially in the 
Aubry-Andre model as well, is essentially a wave phe- 
nomenon and dephasing can be identified as a general 
mechanism of wave packet spreading. As soon as non- 
decaying random phase fluctuations occur localization is 
destroyed with the onset of normal diffusion. Random 
phase fluctuations can also occur due to deterministic 
chaos and nonintegrability of nonlinear wave equations. 
Then the additional dependence of the effective diffusion 
constant on the wave packet density results in a slower 
sub-diffusion as considered in [24] . To maximize the wave 
packet extent at a given time there is an optimal rate and 
strength of dephasing resulting from a competition be- 
tween quickly decohering the initial wave packet and not 
too much slowing down norm radiation to the exterior 
chain. 
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